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Abstract 

A method of cut-off regularization is proposed to evaluate vacuum correc- 
tions in nuclear matter in the framework of the Hartree approximation. Bulk 
properties of nuclear matter calculated by this method are a good agreement 
with results analyzed by empirical values. The vacuum effect is quantitatively 
evaluated through a cut-off parameter and its role for saturation property and 
compressional properties is clarified. 
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In a relativistic approach for the study of nuclear matter vacuum effects 
are very important corrections which can not be taken into account in any 
non-relativistic formalism. About twenty years ago, using a simple a-io model 
(Walecka model), the mean field theory (MFT) clarified the saturation mech- 
anism of nuclear matter relativistically[l] and in succession vacuum effects to 
the saturation property were evaluated by a method of dimensional regulariza- 
tion in the framework of the Hartrcc approximation (RHA)[2]. It was shown 
that vacuum corrections made an effective nucleon mass larger and the incom- 
pressibility of nuclear matter smaller. A small incompressibility is desirable for 
experimental findings of today[3,4]. Recently it has been reported that there 
is a strong correlation between an increase of effective mass and a decrease in 
incompressibility [5] . 

Nearly ten years after a proposal of the Walecka model, corrections due to 
vacuum polarization in quantum fluctuation around the mean field (the Hartree 
field) were evaluated [6]. The results were undesirable. Corrections were larger 
by far than the magnitude of mean field and furthermore brought forth insta- 
ble ghost poles (the Landou ghost) in meson propagators which made nuclear 
matter unstable. One of ad hoc but powerful recipes to escape from these dis- 
asters was to introduce some form factors at each vertex[7]. There was another 
idea that the form factors should be derived from vertex corrections [8]. These 
recipes are grounded on existence of internal structure of hadron. We are afraid 
that vacuum corrections evaluated by the method of dimensional regularization 
may be overestimated because the size of hadron is not considered in such a 
regularization as the one used well in the elementary particle physics. 

When we use a conventional form factor in a calculation of polarization in- 
sertion of vector meson, however, we need a safety device to assure the baryon 
current conservation. In this report, then, we show another recipe to esti- 
mate vacuum corrections in the framework of Hartree approximation in simple 
Walecka model. This is the first request to a method of cut-off regularization 
stated in the following discussion. 

A nucleon propagator G(k) in the relativistic Hartree approximation has the 
following standard form 

G(k)=G F (k) + G D (k) 

= -^--+( - + ^ - - (i) 



E l = + M* 2 , (2) 

M* = M + S s , (3) 

fc* = (fc, fc 4 + £4) = (fc, iko + iSo), (4) 

where the subscripts "F " and "£)" of G(k) denote the Feynman part and the 
density part, respectively, M and fc^ denote the physical nucleon mass and Fermi 
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momentum, respectively. In the simple Walecka model the nucleon self-energies 
are given as follows, 

S S = S., D + E SF = ! A(|L) 2 | 0- t ,{ Gl ,( q ) + G F (i)}- (5) 

where pb denotes the baryon density and m s , m v , g s and g v denote the cr-meson 
mass, the w-meson mass, the a-nucleon coupling and the w-nucleon coupling, 
respectively, and A denotes the degeneracy, A = 2 for nuclear matter and A 
= 1 for neutron matter. The Feynman part of self-energy Y>sf is a divergent 
integral in the 4-dimensional momentum space while T,sd and So have finite 
values. In the Hartree approximation there appears another divergent integral 
in the Feynman part of baryon energy density defined as follows, 

sb=£bd+£bf = ^ j k 2 Vk 2 + M* 2 dk - A j ^—^t Y [ lA G F (fc)] 74 . (7) 

We discuss how we introduce a cut-off parameter into the two divergent 
integrals in the formalism of the Hartree approximation. We start discussion 
with the nucleon effective mass M*. The effective mass is given by using the 
scalar self-energy with a cut-off parameter A, 

M* = Mo + X SD (M*, PB ) + S S f(M*, A) 

= M + E S d{M*, Pb ) + S SF (M*, A) - E S f(M, A), (8) 

M = M-S SF (M,A), (9) 

where Mo denotes the bare mass of nucleon. The last term in eq.(8) is the 
self-energy at zero density and is introduced for the effective mass M* to be the 
physical mass M at zero density. We require the bare nucleon mass M to have 
cut-off dependence so that the physical nucleon mass M dose not depend on 
any cut-off parameter A. We note that the bare nucleon mass M (A) becomes 
the physical mass as A — > 0. 

To make progress we write the energy density as follows, 

e =± 7 r 2 a W >| + A — (M*-M) 2 +sb, (10) 
2 Ztt^ as 

where as and ay arc defined, 

= I'. «^( A ) 2 (11) 

As the second request we require Hugenholtz-von Hove theorem[9], 

(s + P)/p B =E N (k F )=E F , (12) 
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where P denotes pressure and En^f) denotes a nucleon energy at Fermi mo- 
mentum. The condition to satisfy this identity is given by 



de 
dM* 







M* =M 



2 

7r as 



de B 
dM* 



S s = -7r as 



det 



dM*' 



(13) 



The relation Eq.(13) between Eg and Eb is exactly satisfied both in the 
density dependent part and in the Feynman part, respectively, in RHA[2]. In 
the method of cut-off regularization, however, the same relation is not satisfied 
in Feynman part although it is satisfied in the density part. Then we require 
that T,sf is combined with ebf by Eq. (13) and obtain the expression for ebf 
as follows, 



£bf 



=-— [ 



A 1 
44^2 



A 4 ln{ 



M 

A 2 + M* 2 
A 2 + M 2 
A 2 + M 2 



H S f{x,A)-'Ssf(M,A) 



dx 



A 2 _i_ ]i/f *2 

} - A 2 (M* 2 - M 2 ) + M* 4 ln{ + 4 

A 2 + M 2 



} 



^ 4 W^ L }+4(M*-M)M{A 2 -M 2 ln{^^}}] 



where 



Y>sf(M*,A) = -a s M* A 2 - M* 2 ln{ 



^ A 2 



M* 



M* 



(14) 



(15) 



and the lower limit M in the integral is chosen so that Ebf vanishes at zero 
density. The scalar self-energy and the baryon energy density in the density 
dependent part are given by 



^SD = ~a s M* 



k F E* F - M* 2 ln{- 



El 



M* 



£bd 



A 
8^ 



2k F E* F 3 - M* 2 k F E* F - M* 4 ln{ 



k F + E* F 
M* 



(16) 
(17) 



respectively, where E* F — \Jk 2 F + M* 2 . 

We make use of the saturation condition to determine the cut-off parameter 
at the normal density because the cut-off parameter is independent of density. 
From Hugenholtz-von Hove theorem the saturation condition is expressed by 



P = p B [E* f + TT 2 a v p B - e] 
e = M- 15.75 [MeV], 



o. 



(18) 
(19) 



at the normal density. Then we have a following relation between M* and ay 
at the normal density. 



E F +ir 2 a v p B = M - 15.75. 



(20) 
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On the other hand, the coupling strength as is given by 

M* - M 

as = ^ , (21) 

where E s is defined as S s = 0:5 • E s . Putting these as and ay into energy 
density equation, Eq. (10), we have a kind of relation between the effective 
mass M* and the cut-off parameter A which satisfies the saturation condition 
at the normal density. 

In Fig. 1, we show the saturation curves for several sets of M* and A. The 
other bulk properties of nuclear matter can be calculated for each set and are 
summarized in Table. 

Fig. 1 and Table 

First of all we discuss the magnitude of cut-off parameter A which looks very 
small at first glance. We try to calculate T,sf again by introducing the following 
conventional monopole type of form factor into vertex, 

, , A 2 N -M 2 A 2 N -M 2 A 2 
FNN ° M ^Ta^ ■ ( P + q ? + A 2 N ' TTAf' (22) 

Fig. 2 

We note that we use this form factor for T,sf(M, A) at zero density while we 
modify it for S SF (M*, A) at finite density by replacing M with M* in Eq.(22). 
We make a comparison between the new result and S si ?(M*,A), and obtain 
the following relation between two kind of cut-off parameter, 

A^=M 2 + A 2 , (A^ 2 = M* 2 +A 2 ), (23) 

where another cut-off parameter A a dose not participate in Hartrec calculation 
because q 2 = 0. When we take about a half of nucleon mass as the cut-off 
parameter A, for example, we have familiar values for the cut-off A at of the 
monopole type of form factor. 

Next, we are very interested in the saturation curves in Fig. 1. The sat- 
uration curves are softer as an increase of A. Starting from the most stiff 
curve (MFT) obtained without vacuum effect (A = 0), there exist the curves 
with smaller values of incompressibility less than 200 MeV. We can understand 
the saturation mechanism as a stable balance of three force, i.e., the attractive 
a- meson, the repulsive w-meson and another repulsive vacuum effect. Contribu- 
tions of meson to energy density depend strongly on the baryon density whereas 
contributions of vacuum effect depend weakly on the baryon density. So, since 
the energy density is fixed at the normal density, if the vacuum contribution to 
the energy density is larger, the w-meson contribution is smaller. Thus, the in- 
crease of A dulls the density dependence of energy density in the neighbourhood 
of the normal density. This is the reason that the incompressibility K becomes 
small if A increase, as shown in Table. 
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Also wc make a remark on the relation between M* and A. The effective 
mass depends on the attractive cr-nucleon coupling g s and the repulsive vacuum 
parameter A. The former makes M* small and the latter dose M* larger. The 
effective mass continues to increase if < A/M < 0.4, reaches to the maximum 
at A/M - 0.4 and decreases if 0.5 < A/M < 0.56. Then, the vacuum effect for 
the effective mass is largest at A/M <~ 0.4. It is reasonable that the repulsive 
aj-nucleon coupling strength g v shows a minimum value at A/M ~ 0.4 ( the 
maximum effect of vacuum ). It can be observed through the effective mass 
that the parameter A and g v are complement each other as the two repulsive 
effects. 

The effective mass M* decreases as the incompressibility K decreases less 
than 400 MeV. The small M* is rather desirable since the empirical spin-orbit 
splitting in light nuclei supports M* = 0.6M [10]. The skewness K' ( the third 
order derivative of saturation curve at the normal density in Ref. [5]) can also 
be calculated. Using this K', the Coulomb coefficient K c is given by 

based on the scaling model[ll], where g e denotes the proton electric charge. K c 
is the coefficient of leptdermous expansion [11], 

K(A, Z) = K + K.fA- 1 ' 3 + K VS I 2 + K C Z 2 A- 4 / 3 + ■■■, I = 1 - 2^, (25) 

where K s f and K vs are the surface-term coefficient and the volume-symmetry- 
term coefficient, respectively. These coefficients arc determined from the giant 
monopole resonance(GMR) data of many nuclei. In fig. 3, we show the K — K c 
relation together with results analyzed by empirical values in Table 3 of [3] and 
in Table IV of [4] . We have a fine agreement with results analyzed empirically 
if 200 MeV < K < 350 MeV. 

Fig. 3 

The symmetry energy 04 in Table includes the p-meson contribution which 
depends on the square of the ratio of coupling strength to mass of p-meson and 
makes (Z4 increase to about 30 MeV, i.e., without this contribution 04 becomes 
20.6 MeV at K = 300 MeV, in Table. 

The K — K vs relation is shown in Fig. 4. The quantity K vs is given by 

K vs = K sym - L (9^ + 6) , (26) 

in the scaling-model, where L and K sym are the first and second order derivative 
of asymmetry energy, respectively ( see the detailed definitions in Ref. [5] ) . Also 
we have another fine agreement with empirical values[3,4] if 250 MeV < K < 
400 MeV. 

In summary, we proposed the method of cut-off regularization to evalu- 
ate the vacuum corrections in nuclear matter in the framework of the Hartree 
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approximation. We found that this method, RHAC( Relativistic Hartree Ap- 
proximation with Cut-off regularization) , can prepare the values from 200 MeV 
to 546 MeV for nuclear incompressibility in spite of a few adjustable parame- 
ters. So we note that the RHAC method is a very useful phenomenological one 
under the present situation that there is much uncertainty in the experimental 
determination of compressional properties. 

We made quantitative analysis of the vacuum correction by the cut-off pa- 
rameter. The results are summarized as follows. 

(1) An increase of A means an increase of vacuum correction. The parameter 
A can be connected with the cut-off parameter of the conventional monopolc 
type of form factor. 

(2) The vacuum correction gives the repulsive effect both to the effective 
nucleon mass and to the baryon energy density. The saturation property is 
yielded by the interplay among the attractive tr-meson, the repulsive w-meson 
and the repulsive vacuum effect. The repulsive vacuum effect makes the nucleon 
incompressibility small because of its weak dependence on the baryon density. 

(3) The calculated asymmetry energies in Table agree well with the empirical 
values if C 2 p = (g p M /m p ) 2 = 54.71 is used as the p-meson coupling strength[5]. 

(4) The calculated curve on the K — K c plane is a fine agreement with the 
empirical candidates in region 200 MeV < K < 350 MeV and also the curve 
on the K — K vs plane is a good agreement with the empirical candidates in 
region 250 MeV < K < 400 MeV. Therefore, to account for K C ,K VS and a 4 
simultaneously, the RHAC method is valid if 250 MeV < K < 350 MeV. 
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Table and Figure Captions 

Table The numerical results for parameter sets of (M* , A). K, K', K c , 04, 
K vs , K sym and L are shown in MeV. The 04, K vs , K sym and L depend on value 
of C p = g p M/m p where C 2 p = 54.71 [5]. 

Fig.l The fci?-dependence of binding energy. All curves take a minimum value( 
-15.75 MeV ) at the normal Fermi momentum ( 1.35 fm _1 ). The bold solid line, 
the bold dash-dotted line, the bold dashed line, the bold dotted line, the dash- 
dotted line, the dashed line, the dotted line and the solid line are the results 
for K = 546, 500, 450, 400, 350, 300, 250 and 200 MeV in Table , respectively. 
The bold solid line is the result of MFT without the vacuum effect. 
Fig. 2 Diagrammatic representation of the form factor at NNa vertex. The 
solid ( wavy ) lines denote nucleon ( a- meson ) propagators. The form factor is 
shown as the circle at vertex. 

Fig. 3 The K — K c relation. The crosses with error bars are results in ref. [3] 
and the solid squares are the data from the table IV in ref. [4]. 
Fig. 4 The K — K vs relation. The crosses with error bars are results in ref. [3] 
and the solid squares are the data from the table IV in ref. [4]. 
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A/M 


K 


K' 


K c 


a 4 


K 


K 

^^-syrn 


L 


M*/M 


9s 


9v 


0.0000 


546 


226.39 


-8.971 


30.2 


-890.4 


86.3 


100.36 


0.5470 


10.438 


12.899 


0.2827 


500 


160.77 


-8.330 


29.5 


-775.3 


72.0 


95.27 


0.5722 


10.447 


12.498 


0.3884 


450 


101.96 


-7.676 


29.0 


-671.5 


65.6 


91.69 


0.5924 


10.779 


12.165 


0.4892 


400 


54.71 


-7.058 


29.0 


-587.6 


79.8 


92.29 


0.5919 


11.903 


12.173 


0.5407 


350 


0.05 


-6.118 


29.7 


-463.3 


127.9 


98.50 


0.5650 


13.467 


12.615 


0.5559 


300 


-78.33 


-4.320 


30.5 


-197.4 


186.0 


105.02 


0.5398 


14.535 


13.011 


0.5609 


250 


-177.10 


-1.242 


31.1 


288.0 


246.4 


110.74 


0.5202 


15.262 


13.308 


0.5624 


200 


-293.20 


3.972 


31.6 


1142.4 


308.9 


115.86 


0.5044 


15.801 


13.542 



Table 
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